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Abstract 

We demonstrate that the presence of a velocity-dependent term in 
the phenomenological optical potential simulates a source of nonlocal- 
ity. This is achieved by showing that, in the interior of the nucleus, 
the nonlocal wave functions are different from the corresponding local 
ones obtained in the absence of the velocity-dependent term in accor- 
dance with the Perey effect. It is also shown that the enhancement or 
suppression of the nonlocal wave function is energy as well as angular 
momentum dependent. The latter is in line with the results of previ- 
ous works that introduced parity dependent terms in the conventional 
optical potential. 
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1 Introduction 



It is a well known fact that the nonrelativistic nucleon-nucleus optical po- 
tential is nonlocal and energy dependent [H [2]. Nevertheless, a common 
approach to studying the elastic NA scattering is the use of phenomenolog- 
ical local optical models, and a variation of the model parameters with the 
incident projectile energy can be taken as a sign of the presence of nonlocal 
effects in the scattering process [3]. For example, a recent work presented 
local and global optical potentials for neutron and proton elastic scattering 
from nuclei that fall in the mass range 24 > A> 209 while nucleon inci- 
dent energies ranged from 1 keV to 200 MeV. The authors obtained excellent 
elastic angular distribution fits using local optical potentials each consisting 
of 20 fit parameters, however, the strength of the central real part showed 
the largest variation with incident energy. Another measure of nonlocality 
is given by the size of the Percy damping factor, as will be discussed in the 
present study for the case of the presence of a gradient term in the optical 
potential. Such a gradient term, also denoted as velocity dependent term, is 
usually absent in conventional optical model potentials, but has been intro- 
duced in a previous study |5], and will be the object of further analysis in 
the present investigation. 

The nonrelativistic optical potential contains several sources of nonlocal- 
ities. One such nonlocality is due to the Pauli exclusion principle and can be 
taken into account by antisymmeterizing the wave function as in the Hartree- 
Fock theory [6]. Further, in the Hartree-Fock equation, the nonlocality due 
to exchange effects can be expressed in terms of a spatially variable mass 
[7]. More recently, the importance of accounting for the Pauli principle at 
incident energies greater than 25 MeV has been established [S]. At low ener- 
gies the importance of the Pauli exclusion principle has been investigated in 
collective- model coupled- channel calculations [9] . A second source of nonlo- 
cality is due to the coupling of the inelastic excitations to the nuclear ground 
state during the scattering process [2], that gives rise to a nonlocal dynamic 
polarization potential in the elastic channel. Since this nonlocality is diffi- 
cult to take into account rigorously, it is usually accounted for by numerically 
coupling a few important inelastic channels to the elastic channel [T0| [TTj . 
However, when the energy of the incident nucleon is low enough, the scatter- 
ing process is strongly affected by the individual discrete states of the target 
nucleus. A recent study that employed the multichannel algebraic scattering 
theory (MCAS) to explore the dynamic polarization potential [31 [12] for low 
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incident energies also showed that the resulting optical potential was strongly 
nonlocal and energy dependent. 

Recently, we introduced a novel gradient term in the Schrodinger equation 
describing the NA elastic scattering process. We interpreted this term as a 
change in mass of the incident nucleon when moving inside the nuclear matter 
of the target nucleus [5]. The mass change could be a consequence of the 
nucleon-nucleon interaction being affected by the presence of surrounding 
nucleons, like a surface tension which is present only at the surface of the 
target nucleus. This view is supported by the fact that our model, in Ref. 
[5], resulted in excellent fits to the — C elastic scattering data when 
the coefficient of the derivative term was assumed to be proportional to 
the gradient of the nuclear density, which is most important at the nuclear 
surface. Furthermore, our model reproduced well the large angle scattering 
minima that are usually associated with nonlocalities [131 113 • However, it is 
still not yet clear whether the gradient term also simulates other sources of 
nonlocal effects in the phenomenological optical model like, for examples, the 
Pauli exclusion principle and coupling of the elastic channel to the inelastic 
ones. 

We formulated the mass dependence by introducing a Kisslinger type 
potential into the conventional optical model [15] . This is achieved by writing 
the Schrodinger operator as 

where the resulting potential V is velocity-dependent and has the form: 

V{r,p) = 1/(r)H V-p(r)V 

2mo 

^2 r 

= V{r) + p{r)V^ + Vp{r) -V . (2) 

2mo L ^ 

The velocity-dependent potential V{r, p) plays the role of a nonlocal potential 
V{r,r'). In fact there is an equivalence between a velocity-dependent (or 
momentum-dependent) potential and a nonlocal potential V{r,r') [T]. This 
can be seen by interpreting the gradient term that acts on the wave function 
as the first term of a Taylor series that displaces the wave function from point 
r to a different location. Furthermore, the spatially variable effective mass 



3 



■m*(r) is assumed to have the 



1^(1 -pW) (3) 



m*{r) mo 

where p(r) is an isotropic function of the radial variable r that expresses the 
change of the projectile's mass with distance r from the center of the tar- 
get nucleus [5]. The notion of an effective mass has also been introduced in 
relation to relativistic optical models that give rise to significantly reduced ef- 
fective masses [16]. For example, the relativistic effects on the nucleon mean 
free path in the intermediate energy range were investigated in Ref. [17]. In 
addition, the idea of an effective mass was considered for the specification of 
a microscopic optical model potential Ref. [S]. At this point we stress that 
the gradient term in our model is different from the Darwin term present in 
the nonrelativistic equivalent of the Dirac-based relativistic formulation [18] 
because, firstly, our effective mass is real while the Darwin term is complex. 
Secondly, in the nonrelativistic equivalent of the Dirac formulation, the Dar- 
win term is closely coupled to the spin orbit potential while our gradient 
term is independent of the spin orbit term [5]. Furthermore, unlike the case 
for the relativistic models mentioned above, our gradient term is expressed 
as an effective mass that has a radial dependence, rather than in terms of an 
additional scalar potential that is added to the mass. Finally, the effective 
mass formulation has also been introduced in the field of condensed matter 
physics calculations [T ^ [201 1^ . 

Other phenomenological forms of nonlocalities have also been considered 
such as a parity-dependent term in the optical potential [22l |23] for low 
energy scattering of protons from ^^O. This parity dependence was later 
justified tentatively in terms of the Feshbach channel coupling effect [131 [21] ■ 
A presently much used nonlocality, first introduced by Frahn and Lemmer 
[25] and further developed by Percy and Buck in the 1960ies, is given in the 
form of an integral kernel K{r, rf) that acts on the wave function. Percy 
and Buck also obtained a local equivalent potential in the corresponding 
Schrodinger equation [2B] . 

A non local potential V{r,r') acts on the nonlocal wave function ipNL in 
the form J V{r,r')ipNL{r')dr' , while a local potential V{r) acts as a factor 
V{r) ipLir), where is a local wave function. Like the local wave func- 

tion, ipNL is a function of one radial coordinate only. In the nuclear interior, 
a local potential leads to local wave functions which are different from the 
corresponding nonlocal ones obtained using a nonlocal potential, even though 
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both forms may lead to the same asymptotic phase shift. This is known as 
the Perey effect [B], and may be expressed quantitatively as 

ipNLir) = F{r)Mr), (4) 

where ^/'^(r) and ipNLi"^) are the local and nonlocal wave functions respec- 
tively, while the Perey damping factor is defined as F(r)^ [27]. In the nuclear 
interior, an attractive nonlocal potential results in a reduction in the nonlocal 
wave function (-F(r) < 1), while for a repulsive potential the situation is re- 
versed leading to an enhancement of il^Nhij), that is, F{r) > 1. As discussed 
by Austern [6] , nonlocality arises when one singles out from the overall wave 
function that part which describes the motion in the elastic channel. Outside 
the nucleus ipNiif) is the full physical wave function of the system. However, 
inside the nucleus part of the wave function of the system is hidden in the 
invisible (inelastic) channels, which results in a suppression of the nonlocal 
elastic wave function compared to the local one. Conversely, for a repulsive 
nonlocal potential, the particles of the target nucleus are driven away from 
the incident particle. However, conservation of the flux demands that the 
incident particle carries more flux in order to compensate for the part of the 
channel space which has been suppressed [6]. A more recent work employed 
nonlocal potentials in three-body direct nuclear reaction calculations and an 
important nonlocality effect in transfer reactions was reported [28] . 

In this work our aim is to show that the presence of a velocity-dependent 
part in the conventional, phenomenological optical potential simulates a 
source of nonlocality. This is achieved by showing that, in the interior of 
the nucleus, the nonlocal wave functions are different from the correspond- 
ing local ones obtained in the absence of a gradient term (p(r) = 0) in 
accordance with the Perey effect. We shall also demonstrate that this effect 
is angular momentum dependent. 

2 Nonlocality due to the gradient term 

In Ref . [m , we considered a novel velocity-dependent optical potential of the 
form given in Eq. ([2]) based on the change of mass of the incident nucleon, and 
obtained excellent fits for the elastic N —^"^C angular distribution data in the 
low energy range of 12 — 20 MeV. Most notably the model reproduced well 
the prominent large-angle, backscattering minima which depend sensitively 
on the incident energy, and which are a sign of the presence of nonlocalities 
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The local potential consisted of central and spin orbit potentials 
and was taken to have the form: 

V{r) = -Vofir,xo)+Ata^W^^^^+20 [^^\vso+^Wso)l^^^^a■I, 

where Xq stands for [vq, ao) and so on for the rest of the terms. In addition, 
the function f[r,rj,aj) has the following Woods-Saxon form: 

'''' " l + exp[(r-r,AV3)/a^.]' 

with A being the mass number of the target nucleus. In order to reduce 
the number of fit parameters (12) we used the same set of parameters for 
the real and imaginary parts of the spin-orbit term. In what follows the 
nonlocal wave function ipNii''') is the solution of the Schrodinger equation 
corresponding to the velocity- dependent potential V{r,p) given in Eq. (|2]). 
The potential parameters are adjusted to obtain a best fit to the experimen- 
tal N —^"^ C elastic scattering angular distribution. The local wave function 
ipiif), however, is the solution of the Schrodinger equation corresponding 
to the equivalent local potential V{r) whose parameters are adjusted to (i) 
fit experimental data and (ii) produce the same phase shift values as those 
obtained using the nonlocal potential V{r,p). In the local equivalent poten- 
tial case, the effect of nonlocalities in the scattering process are manifested 
as a variation of the potential parameters with energy. This procedure is 
not unique, but is guided by the conventionally used method to fit angular 
distributions. The local wave function obtained in this way is not equal to 
the trivially equivalent local wave function, which is obtained by starting 
the fit in the presence of a velocity dependent term, and then mathemati- 
cally transforming the d il^NL{f)/dr away by a renormalization of the wave 
function {I - p{r)Y/'^^NLir) [29]. 

In Sec. [3l we demonstrate that the nonlocal wave function ipNiif) is in- 
deed different form the local wave function '^/'^(r) in the nuclear interior, 
which supports our hypothesis that the gradient term introduced into the 
conventional phenomenological optical potential simulates a source of non- 
locality. In addition, we find that the enhancement or suppression of the 
nonlocal wave function relative to the local one to be angular momentum 
dependent. This is in line with earlier works that modeled the NA elas- 
tic scattering by introducing a parity-dependent term (—1)' into the optical 
potential 
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3 Results and Discussion 



As described above, the parameters of the equivalent local potential V{r) 
given in Eq. ([5]) were adjusted until (i) a best fit for the elastic — C 
angular distribution data was obtained and (ii) the resulting phase shifts 
were almost identical to the ones obtained corresponding to the nonlocal fits 
in Ref. |5]. The fits for the 12, 14 and 16 MeV incident nucleon energies are 
shown in Fig. [H where the experimental data are taken from Ref. [30] . The 
equivalent local potential fit parameters as a function of the incident energy 
and orbital angular momentum are shown in Table I. 
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20 + 15 i 


0.91 


0.09 


20 


41.0 


1.61 


0.51 


8.0 


1.36 


0.40 


20 + 5 i 


0.80 


0.15 



Table 1: Conventional optical model fit parameters for the — C elastic 
scattering angular distribution. The potential terms are given by Eq. (|5]). 



Clearly, the overall behavior of the angular distributions has been repro- 
duced. However, the fine details of the differential cross sections are less 
well described compared to the corresponding fits obtained using the phe- 
nomenological velocity-dependent optical potential [5]. The large angle back 
scattering minima (associated with nonlocalities) at 18 and 20 MeV were 
poorly reproduced by the local optical potential and hence are not shown. 

Evidently, the fit parameters are energy sensitive, which indicates the 
presence of sources of nonlocalities. Most notable among them being the 
variation in the strength of the central potential Vq. The corresponding non- 
local (velocity-dependent) optical potential fit parameters (given in Tables 
I and II of Ref. |5]) are less energy dependent indicating that the intro- 
duced gradient term only accounts for part of the nonlocalities present in the 
optical model. In a furure work, we intend to investigate other sources of 
nonlocalities in the phenomenological optical potential. 
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In this work we have also investigated the suppression and enhancement 
effects of the nonlocal wave function ipNL{f) as a function of the incident 
nucleon energy in addition to the orbital angular momentum. The Percy 
factor F{r) defined by Eq. (jl]) has been determined by comparing the maxima 
of the real parts of iI)nl{t) and ipiif) in the nuclear interior. Since the 
gradient term simulating the nonlocality is purely real, the imaginary part of 
ipNiif) is almost identical to the corresponding local one. Further, the local 
wave function has been normalized such that it coincides with the nonlocal 
one in the asymptotic region. The damping (or enhancement) factor F{r)^ 
as a function of the incident energy as well as the orbital angular momentum 
is shown in Table [2J 
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1.00 
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1.00 


1.00 


16 


1.134 


0.882 


1.171 


0.918 


1.00 


1.00 


18 


1.355 


0.636 


1.227 


2.001 


1.00 


1.00 


20 


0.928 


0.700 


0.988 


1.522 


0.895 


1.00 



Table 2: Percy factor F(r)^ as defined in Eq. (j4]). The values are calculated 
at the peaks of the nonlocal and local wave functions in the nuclear interior. 

By inspecting Table [2] it is evident that the suppression and enhancement 
effects are both energy and angular momentum dependent. Apart from the 
20 MeV incident energy case, where the angular distribution fits are poorer 
than those obtained at lower energies, the s-wave nonlocal wave function is 
enhanced. For p-wave scattering a suppression is noted for all incident en- 
ergies. Corresponding to / = 2, however, we have a mixed situation where 
the nonlocal wave function is suppressed for low energies and enhanced cor- 
responding to the higher 16 and 18 MeV ones. However, the 20 MeV case 
stands out again and results in a small suppression. For the / = 3 case, we 
have an enhancement at all energies apart from the 16 MeV incident energy. 
For / > 4, the nonlocal wave function is too small in the interior region of the 
nucleus to feel the effect of the gradient term, which simulates the nonlocal- 
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ity, hence F{rY reduces to close to unity. As an illustration, in Figs. [2] and [3] 
we show the real parts of the local 'ipiij) and nonlocal il^Nhij) wave functions 
corresponding to different angular momentum quantum numbers all at the 
same 14 MeV incident neutron energy. As shown in Fig. [2] (a), for the s-wave 
case the nonlocal wave function (solid line) is clearly enhanced in the nuclear 
interior compared to the local one (dotted line). Further, Figs. |2] (b) and [2] 
(c) corresponding to / = 1 and / = 2 respectively show a suppression rather 
than an enhancement. For / = 3, the maximum of ipNiif) coincides with the 
peak of the gradient term (dash-dotted line), which explains the distortion 
of ipNiif) inside the nucleus as shown in Fig. [3] (a). Finally, for the / = 4 
case, the nonlocal wave function in the nuclear interior is too small to be 
affected by the gradient term that simulates the nonlocal effect and acts in 
the internal nuclear region. Therefore, the two wave functions coincide in the 
interior region as can be seen in Fig. [3] (b). The fact that ^/'^(r) and ipNL{f) 
coincide in the external region (for all angular momenta) is guaranteed by 
the adopted normalization method explained above. In addition, since the 
nonlocality is simulated in terms of a purely real term only the real part of 
the nonlocal wave function is affected, the imaginary part is almost identical 
with the corresponding local one and hence is not shown in the figures. 

4 Trivially local equivalent potential 

The Schrodinger equation corresponding to a velocity-dependent potential 
V{r,p) can be transformed into an equation for an equivalent local but 
energy-dependent one, U{r,E), for a trivially equivalent local wave function 
Xl('") through the following transformation on the nonlocal wave function 



which, as derived in [31], leads to the following energy-dependent term in the 
equivalent local potential; 




Xiir) 



(7) 




1 -p(r)' 



(8) 
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where, as defined in Ref. [5], p is given by 

p(r) = poa, ^/^^{ i + exp[(r-r,AV3)/a,] } ' 

In this case the values of the fitting parameters for the nonlocal potential 
V{r,p), given by Eq. ([2]), and the local but energy dependent one U{r,E) 
are identical. Consequently, ipNiir) and Xiir) share the same values of the 
scattering phase shifts. The ratio of the nonlocal wave function ipNiif) to 
the local equivalent one is, by Eq. ([7]), given by the expression 

F(r) = ^ ^ (10) 



This ratio is energy dependent, but has no variation with the orbital angular 
momentum. The energy dependence of F{r)^ arises from the variation of po, 
Tp and ttp with the incident energy. Table |3] shows the maximum values of 
F{rY calculated as a function of the incident energy. The enhancement or 
suppression effects of the nonlocal wave function depend on the value of p(r) 
as given by Eq. ([9]). Table |3] shows that ipNiif) is enhanced for all incident 
energies since < p(r) < 1.0. 



Ei,b (MeV) F{r) 



12 1.82 

14 1.73 

16 1.57 

18 2.25 

20 2.00 



Table 3: Perey factor F{rY as defined in Eq. ( !T0|) . The table displays the 
maximum values of F{rY as a function of the incident energy. 



5 Conclusions 

In a recent work we introduced a velocity-dependent term into the phe- 
nomenological, conventional optical potential. This resulted in excellent fits 
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for the N —^'^C elastic angular distributions even at large back scattering an- 
gles [S], which is a region known to be associated with nonlocalities [131 E]. 
In the present study it has been shown that the velocity-dependent term 
simulates a source of nonlocality by showing that inside the target nucleus 
the nonlocal wave function ipNii'^) is enhanced or suppressed compared to 
the corresponding local one ipiii^) obtained in the absence of the velocity- 
dependent term. This is in accordance with the Perey effect In Table 
[1] we show the fit parameters for the equivalent local potential V{r) given 
in Eq. ([5]) for the N —^"^ C elastic scattering process. Clearly, the parame- 
ters are energy sensitive, and the strength of the local central potential Vq 
shows the largest variation with incident energy. This indicates the presence 
of nonlocal effects in the scattering process. In this work and in Ref. [5] we 
proposed that the gradient term simulates a nonlocality, which is due to a 
change in mass of the incident nucleon when moving inside the nuclear mat- 
ter of the target nucleus. The corresponding fit parameters for the nonlocal, 
velocity-dependent potential V{r,p) (given in Tables I and II of reference [5]) 
show less variation with incident energy. This clearly indicates that there are 
other sources of nonlocalities like, for example, Pauli exchange and knock out 
processes in addition to the detailed structure of the target nucleus that still 
need to be accounted for. The local potential angular distribution fits are 
given in Figure [H which show that the details of the angular distributions 
are less well described compared to the corresponding ones obtained in Ref. 

We have also shown that the enhancement or suppression effects are an- 
gular momentum dependent. In Table [2] we show the Perey factor F(r)^ as 
a function of the incident energy as well as the orbital angular momentum 
quantum number /. The angular momentum dependence of F(r)^ is sugges- 
tive of the inclusion of a parity dependent term ( — 1)' in the conventional 
optical potential [131 l2l] • Since the gradient term simulating the nonlocal- 
ity is purely real only the real part of i/^Niif') is affected and the imaginary 
part is almost identical with the corresponding one for the local wave func- 
tion. Figs, m and [3] show the real parts of local (dotted) and nonlocal (solid) 
wave functions for different angular momentum quantum numbers / all at 
the same incident energy of 14 MeV. Clearly, F{rY is a function of orbital 
angular momentum. 

Finally, we considered two methods to determine the local wave func- 
tion ^/'^(r) corresponding to the nonlocal one, ipNii"^), that was obtained in 
Ref. [3]. The first corresponds to obtaining the best fit parameters for an 
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equivalent local potential such that ipL{r) and 'ipNiij') are phase equivalent. 
The other obtains a trivially equivalent local wave function, which starts 
by fitting the data in the presence of the velocity-dependent term and then 
mathematically transforming the d ipNL{'r)/dr away by a renormalization of 
the wave function as given by Eq. ([7]). In this case the nonlocal ipMLir) and 
trivially local Xl(^) wave functions share the same fit parameters and, by 
Eq. flTUl) . F(r)^ has energy dependence only as shown in Table [31 Since the 
values of p(r), which is given by Eq. (Q, range from zero to less than unity 
'4'NL{f) is enhanced for all energies. 

In conclusion, the non local nature of the derivative term is justified 
here in terms of the properties of the corresponding Perey Damping Factor, 
which differs from unity by 20 to 30%. Furthermore, in future works we shall 
test the applicability of the introduced velocity-dependent potential to other 
systems like p — O and n — Ca, in order to determine whetherthe surface 
term for p(r) as defined in Eq. (|9]) is still valid for heavier nuclei, or whether 
it should be replaced by a volume term that is proportional to the nuclear 
density itself, rather than its gradient. 
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Figure 1: The conventional optical potential fits for N — C elastic scattering 
at 12 (a), 14 (b) and 16 (c) all in units of MeV. The model parameters are 
given in Tabled! The data is obtained from reference |30j . 
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Figure 2: The real parts of the local (dotted) and nonlocal (solid) wave 
functions corresponding to (a) / = 0, (b) / = 1 and (c) / = 2 angular momen- 
tum quantum numbers. The local and nonlocal wave functions correspond 
to the local and nonlocal potentials given by Eqs. ([2]) and ([5]). All figures 
correspond to ii^ = 14 MeV incident neutron energy. 
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Figure 3: The real parts of the local (dotted) and nonlocal (solid) wave 
functions corresponding to the following angular momentum quantum num- 
bers: (a) / = 3, and (b) 1 = 4. In (a) the dash-dotted peak corresponds to 
p(r) defined in Eq. (jS]). The peaks of p(r) and the nonlocal wave function 
almost coincide, which explains the distorted peak of the latter. The local 
and nonlocal wave functions correspond to the local and nonlocal potentials 
given by Eqs. ([2]) and (jS]). Both figures correspond to £" = 14 MeV incident 
neutron energy. 
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